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Abstract 



The present study is a continuation of a previous one on " hyperel- 
liptic" axisymmetric equilibria started in [Tasso and Throumoulopou- 
los, Phys. Plasmas 5, 2378 (1998)]. Specifically, some equilibria with 
incompressible flow nonaligned with the magnetic field and restricted 
by appropriate side conditions like "isothermal" magnetic surfaces, 
" isodynamicity" or P + B 2 /2 constant on magnetic surfaces are found 
to be reducible to elliptic integrals. The third class recovers recent 
equilibria found in [Schief, Phys. Plasmas 10, 2677 (2003)]. In contrast 
to field aligned flows, all solutions found here have nonzero toroidal 
magnetic field on and elliptic surfaces near the magnetic axis. 

PACS: 52.30.Bt, 47.65. +a, 02.30. Jr 
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1 Introduction and basic equation 



A generalized Grad- Shaft anov equation has been derived in Ref. [T] (Eq. (22) 
therein) to describe axisymmetric magnetohydrodynamic equilibria with in- 
compressible flows. This equation consisting the starting point of the present 
investigation is given by 

(1 - Af 2 )A> - i(M 2 )'|W| 2 + \(Y^p)' 

+R (PsW - yztw^ + T ( T^ ) = (1) 

along with a Bernoulli relation for the pressure, 

P = PsW - + -f-], (2) 
z p 

where Ps{ip) is part of the pressure which depends on ip only, ip being the 
poloidal magnetic flux function. The elliptic operator A* is defined by A* = 
i? 2 V • (V/i? 2 ), M 2 = (F'(il))) 2 /p where F(ifj) is the poloidal stream function 
and p(ip) is the mass density, $(V>) is the electrostatic potential, Q/(pR) is 
the toroidal velocity component and X(if>) = 1(1 — M 2 ) +R 2 F'& where I/R 
is the toroidal magnetic field. R, <f), z are the usual cylindrical coordinates 
with z corresponding to the axis of symmetry. As stated in Ref. PQ the 
surface quantities F(ip), Q>(tp), X(tp), p(ip) and Psi.^) are free functions. For 
each choice of this set of five functions, Eq.(l) is fully determined and can 
be solved whence the boundary condition for ip is given. 

For our further investigation it is convenient to simplify Eq.(l) by intro- 
ducing the following transformation 



«(V0= /"^[l - M 2 (g)] 1/2 dg, (3) 
Jo 



which reduces (1) to 



1/ x^ 



2 1 — M 2 ' 



+r 2 (p(v) x *' f ' y i R \ p ^' )2 y-n u) 
+R {Ps{u) ~ T^W 2 * + T ( T^m^ ) " °' (4) 

where the primes indicate now the derivatives with respect to u but F' = 
dF/dip and $' = d<&/dip. (See previous work e.g. Ref. j2j). Note that no 
quadratic term in |Vw| 2 appears anymore in Eq.(4). 
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The paper is organized as follows: section 2 addresses the question of 
the side conditions while in section 3 the shape of the magnetic surfaces is 
determined near magnetic axis. The conclusions are in section 4. 



2 Side conditions 

Instead of specifying the free functions mentioned above to determine Eq.(l), 
it may be of physical or mathematical importance to introduce side conditions 
on some physical quantities like the total pressure, the magnitude of the 
magnetic field or combinations of them. It is indeed plausible to assume 
isothermal magnetic surfaces in hot plasmas PQ because of the huge parallel 
heat conductivity or to try to eliminate neoclassical effects j2] through an 
isodynamic condition. Such side conditions lead, in general, to an additional 
relation between (Vw) 2 , u and R as already accomplished in section 4 of Ref. 
PP or in Ref. |2j. It turns out that, due to the assumed incompressibility of 
the flow, Eq.(4) as well as the side conditions considered here have quartic 
R dependence on the right hand side, which together with Eq.(4) can be 
expressed as follows 



A*u = -f(u) -R 2 g{u) -R 4 h(u), (5) 

|Vu| 2 = 2[i(u) + R 2 j{u) + R A k{u)}, (6) 

where 

X 2 

/(«) = 2(1 _ M iy (7) 
X&F' 

» («> = - mp) • (8) 

and the other coefficients depend upon the specific side condition chosen. 
Let us consider first the case of isothermal magnetic surfaces already treated 
in Ref. jT] in the variable ip. Note that our present variable u defined in (3) 
is a special relabeling of the variable ip. 
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2.1 Isothermal magnetic surfaces 



Setting the plasma pressure P as a function of u and using Ref. [I] to 
calculate the coefficients entering Eq.(6), one obtains 



= -W^Wy (10) 

m = a - - (T^i. (») 
= 2<rb?5)^^>- < 12 > 

To solve equations (5) and (6) simultaneously we use the method of section 
4 of Ref. which boils down to an ordinary differential equation on each 
magnetic surface 



dz 
dx 



P 
Q 



±\[(g + f)x + Uh + k')x 2 -d} 



|2(i + jx + kx 2 ) -j(g + j')x + \{h + k')x 2 - d 



1 1/2- 



(13) 

where x stays for R 2 and five compatibility conditions for seven free functions 
including the five free functions of Eq. (1). There should be no problem, in 
general, to satisfy those compatibility conditions. The solutions of (13) are, 
in general, hyperelliptic integrals jl], which are not related to known special 
functions unless they can be reduced to elliptic integrals. This occurs, in 
particular, for field aligned flows ($' = 0) considered in Ref. [TJ The 
purpose of this contribution is to find other cases of elliptic reduction with 
nonaligned flows ($' ^ 0). The easiest case of that kind is to annihilate the 
coefficient of the largest powers of x in the denominator of Eq.(13) 

h + k! = (14) 

with h and k both different from zero. 

After introducing J = rzrp an d using (9) and (12), Eq. (14) leads to 

J ~ M 2 (1-M 2 ) J = °> (15) 
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whose general solution is 

M 2 

J = C T^W' < 16 > 

with C > and < M 2 < 1. For C = 0, we recover the case of field aligned 
flows already obtained in Ref. pQ. 

As a byproduct of this investigation misprints have been found in the 
nonumbered equations for k{ip) and g(ip) after Eq.(35) of Ref. They 
should be corrected as follows: 

1 o(&) 2 (F'&) 2 

and 

M 2 P s X&F 1 (M 2 )' 

9 ^ ~ 1-M 2 ^JP> ~ \l -M 2 ) 2J + M 2 (l - M 2 ) ' (18) 

2.2 Isodynamic field or B 2 = function of u 

The setting of B 2 as a function of u (Ref. [T] is used for the calculations) 
leads to the functions i(u), j(u) and k(u) entering the side condition (6) as 

'(«) = -^w-y < 19 > 

J'H = (1-M 2 )[ T + Ir - 1 ^], (20) 

fc( u ) = '— M , (21) 

v ' 2(1 — M 2 ) v 7 

while /(«), and are still given by (7)- (9) since Eq. (5) does not 
change. Now the reduction equation (14) becomes 

whose solution is 

with C > and < M 2 < 1. Again we recover for C = the Palumbo 
solution for field aligned flows PQ 03 Ej. 
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2.3 (P + B 2 /2) = function of u 



Though this side condition is not as relevant to hot plasmas as the previous 
cases, it is of mathematical interest since it induces a "hidden symmetry" 
in the equilibrium equations as discovered in Ref. j^j, which leads to a rich 
class of solutions. After calculating (P + B 2 /2) and setting it a function of 
u, we obtain the coefficients i(u), j(u) and k(u) entering condition (6) as 

X 2 

i( u ) = -, (24) 

V ; 2(1 — M 2 ) y ! 

B 2 X&F' 

j( u) = 2[P-P s + - + JT - J ^], (25) 

k(u) = -\j(u), (26) 

while as before f(u), g{u) and h{u) are still given by (7)-(9) since Eq. (5) 
does not change. It turns out that, in this case, the reduction equation (14) 
is identically satisfied, which is reminiscent of the hidden symmetry of Ref. 
[Hj, so recovering the elliptic solutions found therein. 



3 Behaviour of solutions near magnetic axis 

It is instructive to analyse the properties of all possible solutions of (5) and 
(6) near the magnetic axis. Focussing on (5) and (6) we employ a Cartesian 
system (x, y) centred on magnetic axis, i.e. R = -R(O) + x and z = z(0) + y, 
and expand the u surfaces in x and y around the magnetic axis up to second 
order: 

u — n(0) = ax 2 + bxy + cy 2 + higher orders. (27) 

Also, we expand the flux functions contained in (5) and (6) up to first order 
in u — w(0), i.e. 

t(u)=t(0)+t'(0)(u-u(0)), j(u)=j(0)+j'(0)(u-u(0)) etc, (28) 

and R 2 and R A up to second order in x and y. On the basis of the zeroth, first 
and second order equations thus obtained from (5) and (6) we can determine 
the coefficients a, b and c of u — w(0). 
It turns out that 6 = and 

a 1 1 fc(0)fl 2 (0) 1/2 

c ~ 2 + [ 4 + 7 2 J ' 1 J 
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where 7 = [i'(0)+R(0) 2 j'(0)+R(0) A k'(0)]/(2y/2). For the cases of subsections 
2 and 3 Eq.(29), on account of (21) and (26) implying k(0) < 0, means that 
the magnetic surfaces near the magnetic axis are elliptical with elongation 
being directed toward R. For isothermal magnetic surfaces, however, the 
elongation can be either parallel to R if M 2 > 1/2 or parallel to z if M 2 < 1/2 
via (12). For M 2 = 1/2 the isothermal magnetic surfaces become circular. 
Also, in all three cases the ellipses become circles for field aligned flows, i.e. 
the Palumbo's solution is recovered El 12]- This point was overlooked in 
Ref. [I]. In addition, the strength of the toroidal magnetic field on the 
magnetic axis is proportional to k(0)R 2 (0), which vanishes for field aligned 
flows or Palumbo's solution. 

4 Conclusions 

In conclusion, the reduction of the "hyperelliptic" equilibria with nonaligned 
flows introduced in Ref. [T] has been demonstrated for several side conditions 
permitting the discovery of whole classes of magnetohydrodynamic equilibria. 
All these equilibria have nonzero magnetic field and elliptic magnetic surfaces 
on magnetic axis (see Ref. |B. for a special case). The present reduction is 
based on the annihilation of the fourth and fifth power under the square 
root appearing in Eq. (13). Another possible reduction could be sought by 
finding the conditions under which the fifth order polynomial appearing in 
the denominator of Eq. (13) can be factorized with a double zero leaving a 
cubic polynomial under the square root. This and more details about the 
equilibria demonstrated here is left to future work. 
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